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Abstract. The restricted diffusive pair contact process 2A 3A, 2A (PCPD) 
and the classification of its critical behavior continues to be a challenging open problem 
of non-equilibrium statistical mechanics. Recently Kockelkoren and Chate [Phys. Rev. 
Lett. 90, 125701 (2003)] suggested that the PCPD in one spatial dimension represents 
a genuine universality class of non-equilibrium phase transitions which differs from 
previously known classes. To this end they introduced an efficient lattice model in 
which the number of particles per site is unrestricted. In numerical simulations this 
model displayed clean power laws, indicating ordinary critical behavior associated 
with certain non-trivial critical exponents. In the present work, however, we arrive 
at a different conclusion. Increasing the numerical effort, we find a slow drift of the 
effective exponents which is of the same type as observed in previously studied fermionic 
realizations. Analyzing this drift we discuss the possibility that the asymptotic critical 
behavior of the PCPD may be governed by an ordinary directed percolation fixed 
point. 



PACS numbers: 05.50.-|-q, 05.70.Ln, 64.60.Ht 



The phase transition of the diffusive pair contact process revisited 



2 



1. Introduction 

In the attempt to classify non-equilibrium phase transitions from fluctuating phases 
into absorbing states [HEIISI, the so-called pair contact process with diffusion (PCPD) 
continues to attract considerable attention (for a recent review see Ref. |3j). The 
PCPD is a reaction-diffusion process that describes a single species of diffusing particles 
subjected to binary reactions of the form 

2A~^3A, 2A^0 (1) 

combined with a mechanism which prevents the particle density from diverging. In 
one spatial dimension this process exhibits an unconventional type of critical behavior, 
raising the question whether the PCPD represents a novel universality class of non- 
equilibrium phase transitions 0j. Numerical studies of so-called fermionic lattice 
models, where each site can be occupied by at most one particle, displayed a slow drift of 
the effective critical exponents and so far turned out to be inconclusive. More recently, 
however, Kockelkoren and Chate introduced a particularly efficient one-dimensional 
lattice model, in which the number of particles per site is unrestricted [H]. In contrast 
to previously studied models, where the reactions involve three neighboring sites, the 
binary reactions (0) are carried out at single sites, combined with ordinary nearest- 
neighbor diffusion. Unlike fermionic realizations of the PCPD, the model introduced 
by Kockelkoren and Chate was found to exhibit clean power laws over several decades. 
The estimated values of the corresponding critical exponents jH] 

5 = = 0.200(5), z=^ = 1.70(5), /5 = 0.37(2). (2) 

led the authors to the conclusion that the PCPD represents a novel genuine universality 
class of non-equilibrium phase transitions which differs from all previously known classes. 
Likewise, they investigated various other reaction-diffusion rules. In particular they 
found that the so-called triplet process 3A AA, 3A introduced in ^ shows 
yet another novel type of critical behavior. Combining all these results, Kockelkoren 
and Chate presented a generic classification table for absorbing phase transitions 
that comprises four non-trivial classes, namely, directed percolation (DP), the parity- 
conserving (or voter) universality class, the PCPD, and the triplet process. 

The results of the present paper indicate that this classification scheme may be 
premature and needs to be reexamined. Highly optimized large-scale simulations of the 
model introduced in Ref. jO] reveal that Kockelkoren and Chate, apparently attempting 
to obtain straight lines in a double- logarithmic plot, systematically underestimated the 
critical threshold. Increasing the numerical effort it turns out that their model does not 
show clean power laws, instead it is plagued by the very same type strong deviations 
as observed in previously studied fermionic models. The effective critical exponents are 
found to display a slow drift that probably extends beyond the numerically accessible 
range, most likely approaching an ordinary directed percolation fixed point. 
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2. Definition of the model 

The class of models introduced by Kockelkoren and Chate is best explained in the case of 
the PCPD in one spatial dimension (rule ppl2 in their nomenclature). Consider a one- 
dimensional lattice of L sites with periodic boundary conditions. Each site i is associated 
with an integer number Ui > that represents the local number of particles. The 
dynamics is controlled by a single parameter p G [0, 1] and evolves in two synchronized 
substeps as follows: 

• Diffusion: At first each particle diffuses independently with equal probability to 
one of the nearest neighbors. 

• On-site reactions: In the second substep all sites i, which are occupied by at 
least two particles, > 2, are updated individually. To this end the particles 
at a given site are grouped into = [ni/2\ pairs, where [-J denotes truncation 
to the largest integral value not greater than the argument. Each of these pairs 
independently produces one offspring with probability p'^\ increasing Ui by 1, or 
annihilates otherwise, thus decreasing rij by 2. 

Obviously, when p is small, annihilation dominates, leading to an algebraic decay of 
the density p(t) ~ t^^/^ which characterizes the absorbing phase of the PCPD. On the 
other hand, if p is sufficiently large, particle production dominates and the system 
approaches a stationary active phase. Note that for p < 1 the generation of new 
particles is exponentially suppressed as rii increases. This exponential cutoff, which 
is the counterpart of the fermionic constraint (nj < 1) in previously studied PCPD 
models, prevents the particle density from diverging in the active phase. 

The update algorithm can be optimized as follows. Because of the exponential 
cutoff particle numbers greater than rij = 16 are extremely rare so that we can safely 
restrict their range to < < 32. To speed up the on-site reactions we define for 
each Ui a list of probabilities for the possible outcomes which are then selected by 
using a single random number. Furthermore, instead of implementing the lattice as a 
static array, the actual configuration is stored as a dynamically generated list of the 
coordinates of occupied sites. These two improvements, combined with some general 
technical optimizations, accelerate the update algorithm significantly. 

3. Numerical results 

Let us consider simulations starting with a homogeneous initial state, where each site is 
occupied by a pair of particles. Kockelkoren and Chate studied this case by performing a 
single run on a large system with L = 2^^ ^ 4.2 x 10^ sites simulating up to t = 2 x 10^ 
Monte Carlo updates (see Fig. 1 of Ref. [Hj). Plotting various order parameters in 
a double-logarithmic representation and seeking for asymptotically straight lines they 
estimated the critical point pc = 0.795410(5). In the improved simulations presented 
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Figure 1. Density of particles starting with homogeneous initial conditions as a 
function of time muhiphed by t^''^ for various values of p. The error bars indicate the 
expected standard deviation averaged over the samples. Other order parameters as 
those studied in Ref. lead to similar results. The dashed line represents the slope 
that corresponds to the critical behavior of directed percolation. 

here the temporal range is extended by almost one decade up to t = 1 x 10* Monte 
Carlo updates. Furthermore, using the same system size of L = 2^^ sites as in Ref. jU] 
we average over at least 40 independent realizations of randomness in order to reduce 
the statistical error. The results are shown in Fig. ^ As can be seen, the curves for p 
equal to 0.795405, 0.795410, and 0.795415 bend downwards. Consequently the critical 
point pc = 0.795410(5) given in jH], including its error margin, lies entirely in the inactive 
phase. It seems that Kockelkoren and Chate, attempting to obtain an apparent power- 
law behavior, systematically underestimated the critical point. 

Instead of clean power laws we observe a slow drift of the effective critical exponents 
which manifests itself as a weak curvature of the data in Fig. ^ that extends over at 
least four decades in time. Therefore the model introduced by Kockelkoren and Chate 
is not a good-natured exception in the zoo of suggested PCPD models, rather it shows 
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z 




Ref. 


0.795410(5) 


0.200(5) 


1.70(5) 


0.37(2) 


present work 


0.795417(2) 


< 0.185 


< 1.65 


< 0.34 


DP 




0.1595 


1.5807 


0.2765 



Table 1. Estimates for the critical exponents compared with DP exponents. 

the same phenomenological properties as previously studied fermionic variants and is 
plagued by exactly the same type of slowly decaying corrections to scaling. At the 
present state of knowledge there is no reason to believe that the observed drift suddenly 
stops after 10^ or 10^ time steps, instead it seems to be more likely that the drift 
continues to vanish gradually over many decades beyond the numerically accessible 
temporal range. Therefore, it would be a mistake to search again for straight-looking 
tails of the curves and to specify updated exponents of a possible PCPD universality 
class. However, assuming that the drift eventually vanishes without changing sign, such 
estimates may serve as upper bounds of postulated asymptotic critical exponents. For 
example, analyzing the data shown in Fig. Qone obtains 6 < 0.185. Similar inequalities 
are obtained by finite-size and off-critical simulations, as summarized in Table CI 

4. Discussion 

The numerical results presented in this paper demonstrate that the PCPD model 
introduced in Ref. [Oj, in contrast to earlier claims, does not reach the scaling regime in 
the numerically accessible temporal range, instead it displays a slow drift of the effective 
critical exponents, just in the same way as in previously studied 'fermionic' variants 
of the model with hard-core exclusion. Because of this drift the critical exponents 
estimated by Kockelkoren and Chate turn out to be incorrect and have to be replaced 
by the inequalities given in Table 

As in all PCPD models the observed drift is so slow that any attempt to extrapolate 
the effective exponents to t — oo involves to some extent speculative expectations. 
Strictly speaking it is even not yet clear whether such a limit exists, questioning the 
very existence of a second-order phase transition. However, in my opinion the numerical 
results strongly support earlier suggestions [01 1101 that the one-dimensional PCPD - 
including its realization by Kockelkoren and Chate - may eventually cross over to an 
ordinary directed percolation transition. 

The hypothesis of a DP transition in the PCPD is based on the following argument 
(cf. Ref. in] for details). In the numerically accessible range the PCPD is neither 
reaction- nor diffusion-limited, instead one observes a spatio-temporal coexistence of 
clustered reactions and freely diffusing solitary particles 0. These single particles diffuse 
over large distances in the l-|-l-dimensional geometry with the dynamic exponent z = 2 
and their diffusion constant is well-defined and can be read off easily at any time. 
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On the other hand, all simulations in 1+1 dimensions consistently confirm that the 
critical cluster spreads superdiffusivcly, i.e., z < 2. The assumption of a scale-free 
situation therefore leads to a contradiction unless the diffusion constant flows to zero 
under renormalization group transformations as the critical point is approached, i.e., to 
a situation which is known to belong to the directed percolation universality class. 

In the light of the present numerical results the classification scheme introduced 
by Kockelkorcn and Chate needs to be reexamined carefully by various methods. In 
particular it would be interesting to verify their conclusions concerning the triplet 
process by similar large-scale simulations. 
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